A theoretical and numerical approach to "magic angle" of stone skipping 
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We investigate the condition for the bounce of circular disks which obliquely impacts on fluid 
surface. An experiment [ Clanet, C, Hersen, F. and Bocquet, L., Nature 427, 29 (2004) ] revealed 
that there exists a "magic angle" of 20° between a disk's face and water surface in which condition 
the required speed for bounce is minimized. We perform three-dimensional simulation of the disk- 
water impact by means of the Smoothed Particle Hydrodynamics (SPH). Futhermore, we analyze 
the impact with a model of ordinal differential equation (ODE). Our simulation is in good agreement 
with the experiment. The analysis with the ODE model gives us a theoretical insight for the "magic 
angle" of stone skipping. 

PACS numbers: 45.50.Tn, 47.11. +j, 47.90.+a 



Problem of impacts and ricochets of solid bodies 
against water surface have been received a considerable 
amount of attention 

the problem was of importance in naval engineering con- 
cerning the impacts of canon balls on sea-surface pj- In- 
vestigations then revealed that there exists a maximum 
angle of incidence # m ax for impacts of spheres, above 
which the rebound does not occur 8]. Besides, it was 
empirically found that the # max relates to specific grav- 
ity of sphere a as # max = 18/y/a. This relation was the- 
oretically explained using a simple model of an ordinal 
differential equation (ODE) |8|,|9(. In military engineer- 
ing today, the problem of water impacts may be not as 
important as that of a century ago, however, recently it 
attracts renewed interest under the studies of locomotion 
of basilisk lizards and stone-skip . 

This study is motivated by experimental study of 
stone-skip - bounce of a stone against water surface - 
by C. Clanet et. al. They investigated impacts of a 
circular disk (stone) on water surface and found that an 
angle about <p = 20° between the disk's face and water 
surface would be the "magic angle" which minimizes re- 
quired velocity for bounce. In this paper, we study theo- 
retically and numerically the oblique impact of disks and 
water surface. Our simulation successfully agrees with 
the experiment. Moreover, we apply an ODE model ^3 
to the disk-water impact and obtain an analytical form 
of the required velocity v m [ n and maximum angle # max as 
a function of initial disk conditions. 

To perform a numerical simulation of the disk-water 
impact, we solve the Navier- Stokes equation using the 
technique of Smoothed Particle Hydrodynamics (SPH) 
[HI [3- Fig- Glis the snapshots of our simulation. The 
SPH method is based on Lagrangian description of fluid 
and has an advantage to treat free surface motion. Sev- 
eral representation of the viscous term have been pro- 
posed for this method. In this work, we adopt an arti- 
ficial viscous term [T^j which is simple for computation 
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FIG. 1: Snapshots of the SPH simulation of a disk- water im- 
pact with (j> = 20° and incident angle 6 — 15. Specific gravity 
of the disk a = 1.5 and number of fluid particles N = 12600. 



and sufficiently examined with Couette flow [l6|- In our 
simulation, we neglect surface tension and put the veloc- 
ity of sound of the fluid, at least, 25 times larger than 
the incident velocity of the disk. 

In the following discussion, we analyze the ODE model 
which was originally introduced by Birkoff et al 17]. The 
model is based on the following assumptions; (i) Hydro- 
dynamic pressure p acting from water is proportional to 
(v • n) 2 , where v is the speed of the body and n is the unit 
vector to the surface of the disk, (ii) For the part of the 
surface facing air, there is no hydrodynamic force, (iii) 
During the whole process, deformation of water surface 
is negligible, and the boundary between immersed and 
non-immersed area is simply given as the cross section 
to a horizontal plane at water level. We notice that the 
first assumption is reasonable because Reynolds number 
would be of order 10 5 for typical cases of stone-skip [T^ . 
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FIG. 4: Relation of incident angle and the angle <j>* . The 
SPH simulation is performed with a — 2.0. The solid line is 
obtained numerically seeking the minimum of Eq. Q. 



FIG. 3: A comparison of the stone-skip domains obtained 
from the experiment SPH simulations and our theory, 

(i), The stone-skip domain in the (v, <$) plane for a fixed 6 = 
20°. (ii), The stone-skip domain in the (0, 0) plane for a fixed 
v = 3.5 [m/s]. The boundary of the stone- skip domain under 
the criterion A in each graph are numerically drawn and those 
of B are the plot of Eqs.Q and respectively. 



Note that, in the limit F — > oo, this equation again gives 
cr + <j> = 7r/2. 

A position of the inflection point still remains un- 
known. We treat as a fitting parameter, which should 
be determined so as to agree with experiments. However 
we cannot make a direct comparison between the Eq. (Q) 
and the experimental data (P3) because v m [ n and # max 
are acquired under the criterion A in the experiment. We 
thus fit the Eq. J7J) with the result of our SPH simula- 
tions performed under criterion B and evaluate = 2.6. 
Due to the nature of the criterion B, these analytical ex- 
pressions for v m i n and # m ax should give a lower and upper 
limit of the stone-skip domain respectively. 

Then let us discuss our results. We chose the same 
parameters as that of the experiment: A = 9.1 and 
a = 2.7 unless particularly mentioned. Froude num- 
ber F typically ranges from 4.0 to 200. For the SPH 
simulation, A = 2.5 and angular velocity of the disk 



uj = 65[rounds/s]. Figure 01 shows the domains of stone- 
skip in ((/>, v) and (0, (j>) planes. For the minimum veloc- 
ity ^min, the SPH simulation successfully agrees with the 
experiment, and the theoretical results under criterion A 
and B also show the qualitative agreement. 

The experiment indicates that the stone-skip domain 
shrinks at < 20° in (0, <$) plane. The theoretical curve 
under criterion B does not reproduce this tendency while 
that of criterion A shows the similar behavior. This in- 
consistency is due to the assumption that the disk is fully 
immersed in the water when it reaches to the inflection 
point. However, in the case that the is much smaller 
relative to the tilt angle (f) this is totally incorrect: only 
small part of the disk is immersed during the impact 
process. The SPH simulation also shows the different 
behavior with the experiment under < 20°. We cannot 
present a clear explanation for this discrepancy. As for 
SPH simulations, we mention that, the depth of immer- 
sion of the disk would be of the order of the fluid particle 
size of SPH at very small incident angle. The numerical 
error hence becomes larger for small and for the domain 
< 10° simulation is not attainable. 

The angle <p ~ 20° is a characteristic for both ((/>, v) 
and (0, (/)) planes in the experiment. C. Clanet et. al. 
hence suggested that the angle <\> = 20° is the "magic 
angle" for stone- skip. However Eq. J7J) implies that <\> 
depends on 0. In Fig. we show how the "magic angle" 
(j) m is affected by the incident angle 0. Our theory sug- 
gests m decreases as incident angle increases and SPH 
simulation also shows a decreasing tendency. However, 
the change in (j) m is sufficiently small: (j) m changes only 
about 15% relative to the change of incident angle under 
= 40°. We therefore conclude that the "magic angle" 
still remains around <fi = 20° for the ordinal incident an- 
gle at stone skipping. 
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